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We present branching fraction measurements of charged and neutral B decays to D n~ , D* n~ and 
D**tt~ with a missing mass method, based on a sample of 231 million T(4S) — > BB pairs collected 
by the BABAR detector at the PEP-II e + e~ collider. One of the B mesons is fully reconstructed 
and the other one decays to a reconstructed charged ir and a companion charmed meson identified 
by its recoil mass, inferred by kinematics. Here D** refers to the sum of all the non-strange charm 
meson states with masses in the range 2.2 — 2.8 GeV/c 2 . We measure the branching fractions: 



o 



B(B~ 

B{B~ -> L>* u tt 

— > D TV 

B(B° ^D + tv 
B(B° -> D* + ty 



B(B 



B(B° 



= (4.49 ± 0.21 ± 0.23) x 10~ 3 

= (5.13 ±0.22 ±0.28) x 10~ 3 

= (5.50 ± 0.52 ± 1.04) x 10~ 3 

= (3.03 ±0.23 ±0.23) x 10" 3 

= (2.99 ± 0.23 ± 0.24) x 10" 3 

= (2.34 ±0.65 ±0.88) x 10" 3 



and the ratios: 



B(B~ 
B(B~ - 
B(B° 



■» D*°n-)/B(B- 
D**\-)/B(B~ 
> D* + tv-)/B(B° 



D u n~) = 1.14 ±0.07 ±0.04 
D°-k~) = 1.22 ±0.13 ±0.23 
D + ?v~) = 0.99 ± 0.11 ± 0.08 
D + tv~) =0.77 ±0.22 ±0.29 



B(B° -> D** + n~)/B(B° 
The first uncertainty is statistical and the second is systematic 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Our understanding of hadronic B-meson decays has 
improved considerably during the past few years with 
the development of the Heavy Quark Effective Theory 
(HQET) [US and the Soft Collinear Effective Theory 
(SCET) |3j,|4|. In these models, and in the framework of 
the factorization hypothesis 0, , the amplitude of the 
B — > D(*'ft two-body decay carries information about 
the difference 5 between the strong-interaction phases 
of the two isospin amplitudes A 1 / 2 and A 3 / 2 that con- 
tribute 0, 0- A non-zero value of S provides a mea- 
sure of the departure from the heavy-quark limit and the 
importance of the final-state interactions in the D^tt 
system. With the measurements by the BABAR || and 
BELLE experiments of the color-suppressed B decay 
B° — > Z)(*) 7r providing evidence for a sizeable value of 
<5, an improved measurement of the color-favored decay 
amplitudes (B~ D^ir' and B° D^ )+ tt~) is of 
renewed interest. In addition, the study of B decays into 
D, D* , and D** mesons will allow tests of the spin sym- 
metry [U [H E3 imbedded in HQET and of non- 
factorizable corrections EH that have been assumed to 
be negligible in the case of the excited states D** 

In this paper we present new measurements of the 
branching fractions for the decays B~ — ► D° ir^^D* tt~ , 
D**°ir~ , and B° -> D+ tt~ , D*+ k~, D**+tt- [3, based 
on a missing mass method previously used by BABAR E3 • 
Here D** refers to the sum of all the non-strange charm 
meson states with masses in the range 2.2 — 2.8 GeV/c 2 . 
This analysis uses T(4S) ~ > BB events in which a B + 



or a B° meson, denoted B reco , decays into a hadronic 
final state and is fully reconstructed. The decays of 
the recoiling B into a charged pion and a charmed me- 
son, i.e. B — > tt~X, are studied. The charged pion 
is reconstructed and the mass of the X = D,D*,D** 
is inferred from the kinematics of the two body B de- 
cay. This method, unlike the previous exclusive measure- 
ments does not assume that the T(AS) decays 
into B + and B° with equal rates, nor does it rely on the 
D, D*, or D** decay branching fractions. 

The measurements presented here are based on a sam- 
ple of 231 million BB pairs (210 fb -1 ) recorded at the 
T(4S) resonance with the BAB4R detector at the PEP- 
II asymmetric-energy B factory at SLAC. The BABAR 
detector is described in detail elsewhere [^ij. Charged- 
particle trajectories are measured by a 5-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH), both operating in a 1.5-T solenoidal 
magnetic field. Charged-particle identification is pro- 
vided by the average energy loss (dE/dx) in the track- 
ing devices and by an internally reflecting ring-imaging 
Cherenkov detector. Photons are detected by a CsI(Tl) 
electromagnetic calorimeter. Muons are identified by the 
instrumented magnetic-flux return (IFR). We use Monte 
Carlo (MC]simulations of the BABAR detector based on 
GEANT4 |21| to optimize selection criteria and deter- 
mine selection efficiencies. 

We reconstruct B + and B° decays (B reco ) in the modes 
B+ -> D^ a ir+, D(*)°p+, ^ W V, and B° 
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£)(*)- 



£)(*)- 



the if + 7r~, K + i 



D° 

T-0 



candidates are reconstructed in 



and 7<97r + 7r de- 



cay channels, while D~ candidates are reconstructed 
in the K + n~ir~ and K®tt~ modes, and K® mesons 
are reconstructed to 7r + 7r~. D* candidates are recon- 
structed in the D*~ -> D°ir- and D*° -> ZJ°7r° de- 
cay modes. A 3<r cut is applied on the D meson mass 
rriD (and on the D*-D mass difference Amu*) where 



is the resolution on run (Amc«) and 



is determined from data. A vertex fit is performed on 
D (D*) with the mass constrained to the nominal value 
[2^ . Two nearly independent variables are defined to 
identify the fully reconstructed B candidates kinemati- 
cally. The first one is the beam-energy substituted mass, 
toes = \J (s/2 + p l ■ p B ) 2 /Ef - p|, where p B is the 
-Broco momentum and (Ei,pi) is the four-momentum of 
the initial e + e~ system, both measured in the labora- 
tory frame. The invariant mass of the initial e + e~ sys- 
tem is \/s. The second variable is AE = E* B — y/s/2, 
where E B is the B reco candidate energy in the center- 
of-mass frame. To define the B roco sample (Fig. HJ, we 
require \AE\ < naAE, where the measured resolutions 
(Jae range from 12 to 35 MeV and n = 2 or 3, both 
depending on the -B 10C o mode. The -B roC o candidate mul- 
tiplicity is 1.4 for data as well as for the MC simulation 
sample. For events with more than one candidate, we se- 
lect the B r0 co with the best \ 2 defined with the variables 
mo, Amp*, and AE. The MC simulation shows that 
the recoil variables are reconstructed well within their 
experimental resolution when using this selection. 

The number of -B loco is extracted from the toes spec- 
tra (Fig.Q} in the 5.27- 5.29 GeV/c 2 signal region. The 
toes distribution is fit to the sum of a broad combi- 
natorial background and a narrow signal in the mass 
interval 5.21 — 5.29 GeV/c 2 . The combinatorial back- 
ground is described by an empirical phase-space thresh- 
old function and the signal with a Crystal Ball func- 
tion [24| which is a Gaussian function centered at the 
B meson mass modified to account for photon radiation 
energy-loss. All parameters for the functions describing 
the -B TO co signal and background distributions are deter- 
mined from data. The measured yields of reconstructed 
B+ and B° candidates, N B + = 189474 ± 7487 and 
N B a — 103169 ± 3303, are obtained by subtracting 
the fitted and the peaking (described below) backgrounds 
from the total number of events found in the signal re- 
gion. These i? r0 co numbers serve as the normalization 
of all branching fraction measurements reported in this 
paper. The error is dominated by the systematic uncer- 
tainties due to the fit of the combinatorial background 
and to the determination of the peaking background. We 
assign 2.3% uncertainty to N B + and 1.8% to N B o as a fit 
uncertainty, obtained by varying the lower boundary of 
the fit interval from 5.20 to 5.23 GeV/c 2 . The contam- 
ination of misreconstructed B° events in the B + signal 
(and vice-versa) induces a peaking background near the 
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FIG. 1: m,ES spectra of reconstructed (a) B + and (b) B° 
candidates. The solid curve is the sum of the fitted signal 
and background whereas the dashed curve is the background 
component only. 



B mass. From the MC simulation, the fraction of B° 
events in the reconstructed B + signal sample is found 
to be (3.2 ± 3.2 syst )% and the fraction of B + events in 
the reconstructed B° signal sample (2.8 ± 2.8 syst .)%. A 
100% systematic uncertainty is conservatively assigned 
to these numbers taking into account the possible differ- 
ences in the reconstruction efficiency in data and MC, 
as well as the branching fraction uncertainties for those 
B decay modes contributing to the peaking background. 
The total systematic uncertainties on N B + and N B o arc 
3.9% and 3.2%, respectively. 

In the decay T(4S) — > B ICCO Bxn where Bx-n is the 
recoiling B which decays into 7r~A, the invariant mass 
of the X system is derived from the missing 4-momentum 
px applying energy-momentum conservation: 



PX = PT(4S) ~ Pb ibc . 



Ptv- 



The 4-momentum of the T(4S I ), j>r(4S)i is computed from 
the beam energies and p n - and ps roco are the measured 
4- momenta of the pion and of the reconstructed B Teco , 
respectively. The B TCCO energy is constrained by the beam 
energies. The B -> Dtt~ , B -> D*n-, or B -> D**ir- 
signal yields peak at the D, D*, and D ** masses in the 
missing mass spectrum, respectively. 

The pion candidates, chosen among the tracks that do 
not belong to the B rcco , arc required to have produced at 
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least 12 DCH hits. For the charged B rcco , the pion candi- 
date has the opposite sign to the B [oco . For neutral B reco , 
because of the B a —B° mixing, the corresponding require- 
ment is not applied. Muon tracks are rejected using the 
IFR information, electrons tracks using the energy loss 
in the SVT and the DCH, or the ratio of the candidate's 
EMC energy deposition to its momentum (E/p). Pro- 
tons and kaons are rejected based on informations from 
the DIRC and energy loss in the SVT and the DCH. The 
rejection efficiency is 97% and there is no peaking trend 
in the missing mass distribution from remaining kaons, 
protons, muons, or electrons. The pion reconstruction 
efficiency is determined from the MC simulation and re- 
ported in Table [I] 

The signal yields for the different decay modes are ex- 
tracted from the missing mass spectra. The data distri- 
butions and the bb and the qq (q = c, u, d, s) background 
expectations are shown in Figs. Eta) and|2fb). The shape 
of the background is taken from MC and the normaliza- 
tion is scaled to match the data in the sideband region 
2.8 — 3.2 GeV/c 2 . The error on the background normal- 
ization is 2%. This is determined using the statistical 
errors of MC and data samples. The background sub- 
tracted missing mass distributions are shown in Figs.EIc) 
and [3d). 

The Dtt and D*ir signal yields are extracted by a \ 2 
fit to the background subtracted missing mass distribu- 
tion in the range 1.65 - 2.20 GeV/c 2 . The Dtt and D*ir 
components are each modeled by a sum of two Gaus- 
sian functions, to account for tails in the mass distribu- 
tions. The parameters arc mf ' and a[ >i ) for the D 
and D* resonances, where the index £ = 1,2 corresponds 
to the first and second Gaussian. In the fit, the cen- 
tral values mf and the erf are free parameters, while 
for the D* the variances are constrained by the ratios 
of*/ "? = 0.900 ± 0.015, as determined from MC sim- 
ulation, while the central values differences mf — mf 
arc fixed to 0.1421 GeV/c 2 and to 0.1406 GeV/c 2 for B+ 
and B° , respectively, corres pon ding to the world average 
D and D* mass differences |22j. 

The D** yields are defined as the excess of candidates 
in the missing mass range 2.2 — 2.8 GeV/c 2 , and the 
B — > D**n~ branching fractions refer to the contribu- 
tions of all non-strange charm meson states in the same 
region. The range is chosen in order to maximize the 
acceptance to the four P-wave D** states predicted by 
the theory given the 34 MeV/c 2 mass resolution, deter- 
mined from MC simulation, in the same region. The 
well-known narrow D\ and D\ states are fully con- 
tained in this range, and more than 90% of the broad 
Do and D[ , are covered if measured masses and widths 
|25l l2q are used. The event yields, the efficiencies, and 
the resulting branching fractions are reported in Table [I] 

The uncertainty related to ir reconstruction efficiency 
is due to the MC sample statistics and the systematic 
uncertainty on track reconstruction and particle iden- 



TABLE I: Signal yields, efficiencies and branching fractions 
for B —* Dir~ , B — » D*ir~ and B — » D**tt~ . The first error 
is statistical except for the efficiencies for which it is mainly 
systematic. The second error on the branching fractions is 
systematic. The B — * D**ir~ branching fractions are given 
for the 2.2 — 2.8 GeV/c 2 mass range which in addition to the 
P-wave states may include some yet unknown charm meson 
states. 



Decay mode 



Yield 



Efficiency 



B(10" 



D*°7r- 



677 ± 32 

774 ± 33 0.796±0.007 

829 ± 78 



4.49±0.21±0.23 
5.13±0.22±0.28 
5.50±0.52±1.04 



B" 
B° 
B° 



D+n- 

D*+n- 

D**+tt- 



248 ± 19 

245 ± 19 0.793±0.007 

192 ± 54 



3.03±0.23±0.23 
2.99±0.23±0.24 
2.34±0.65±0.88 



tification algorithms. The uncertainty due to the yield 
extraction is estimated by fitting the MC sample. The 
difference between the MC and the data fitted yields 
is found to be consistent with zero and the statistical 
errors are taken as a systematic error. We evaluate the 
uncertainty on the missing mass resolution in the Dn 
and D*n yield extraction by varying by one standard 
deviation the ratio crP /af while a® and 
are let free. The difference in the yield is taken as 
systematic uncertainty. The uncertainty related to the 
subtraction of the background is determined by varying 
the branching fraction of the different background 
components within the uncertainties of the most recent 
measurements |22l | and taking into account the error on 
the background normalization. Due to the threshold 
shape of some of the background components and the 
fast varying combinatorial background, B — ► D**tx 
branching fractions have larger systematic errors than 
B — ► Dtt and B — ► D*ir branching fractions. The 
summary of these systematic uncertainties is reported 
in Table HD 

Using the measured branching fractions we compute 
the following ratios: 

B(B- -► D*°tt-)/B(B- -> D°tt-) = 1.14 ±0.07 ±0.04, 
B(B~ -> D**°tt-)/B(B- -> D°tt-) = 1.22 ± 0.13 ± 0.23, 
B(B° -► D* + tt-)/B(B° -» D+TT-) = 0.99 ± 0.11 ± 0.08, 
B(B° -> D**+n-)/B{B° -> D+tt") = 0.77 ± 0.22 ± 0.29. 

The first uncertainty is statistical and the second is sys- 
tematic. In addition to the cancellation of many of the 
systematic errors, the ratios are insensitive to the abso- 
lute normalization scale. 

In summary, we have measured the branching frac- 
tions for the decays B~ — > D°ir~ , B~ — ► £>*°7r~, 
B~ -► D**°tt-, B° -► D+Ti-, B° -> D*+ir-, and 
B° — > D** + ir~ , using a missing mass method. This mea- 
surement does not assume that the T(4S) decays into B + 
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FIG. 2: Top: missing mass distributions obtained in the recoil of B + (a) and B° (b). The points with error bars show the data 
and the histograms show the background contributions (bb and qq (q = c,u,d,s)) predicted by the MC simulation. Bottom: 
background-subtracted missing mass spectra for B + (c) and B (d). The curves show the result of the fits to the Dn and D*-n 
components. 



B(B° 



II: Total relative 


systematic 


uncertainties for 


the branching 


; fractions 


B(B~ 


- (D°, 




-> (d+, d*+, d**+)tO 
















Syst. Source 


B- -f D° 


7T — B~ -> D*°7r" 


B - _^ £) **o 7r - 


B u -» 


- B u 




B u -» 73"+^- 


N B 


3.9% 


3.9% 


3.9% 


3.2% 




3.2% 


3.2% 


Efficiency 


0.9% 


0.9% 


0.9% 


0.9% 




0.9% 


0.9% 


Yield extraction 


2.7% 


2.7% 


5.1% 


5.4% 




5.1% 


5.9% 


Missing mass resolution 


0.9% 


0.8% 




1.9% 




1.1% 




Background subtraction 


1.6% 


2.3% 


17.7% 


3.7% 




5.4% 


37.1% 


Total 


5.2% 


5.4% 


18.9% 


7.6% 




8.2% 


37.7% 



and i?° with equal rates, nor does it rely on the D, D*, 
or D** intermediate branching fractions. The results for 
B(B — > Dir~) and B(B — > D*ir~) are compatible with 
previous world averages [12 ■ We have extracted a new 
result for B(B — > Z3**7r~) branching fractions where £>** 
excited states correspond to the yield measured in the 
mass range 2.2 — 2.8 GeV/c 2 . The isospin study 0, Q 
will become competitive with the exclusive measurements 



|l9j if the statistical error is reduced by a factor of 
2. With regard to spin symmetry, the values measured 
forjhe ratios B{B~ D*°ir-)/B(B- -> D°7r-) and 
D*+tt-)/B(B° D+tt-) are close to 1, as pre- 
dicted by different theoretical models [Tol HH HI HI flU . 
and their precision is comparable or better than the cur- 
rent world averages [22^ . 
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